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Abstract 

Reaction paths and classical and quantum trajectories are studied within a generalized Hamilton- 
Jacobi framework, which allows to put on equal footing topology and dynamics in chemical reactiv- 
ity problems. In doing so, we show how high-dimensional problems could be dealt with by means of 
Caratheodory plots or how trajectory-based quantum-classical analyses reveal unexpected discrep- 
ancies. As a working model, we consider the reaction dynamics associated with a Miiller-Brown 
potential energy surface, where we focus on the relationship between reaction paths and trajectories 
as well as on reaction probability calculations from classical and quantum trajectories. 
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I. INTRODUCTION 



One of the most outstanding problems in theoretical Chemical Physics is studying and 
determining the mechanisms that underly chemical reactions [l]. Traditionally, a full the- 
oretical study of chemical reactions involves two steps: (1) determining a potential energy 
surface (PES) and (2) the analysis of the associated dynamics. This has given rise to two 
types of comparative schemes. Within the first scheme, comparisons are established between 
the topological properties of PESs and the dynamics ruled by them. Because most of the 
energy exchange is carried by electrons in the reactants-to-products transition, chemical re- 
actions are usually studied neglecting the (tinny) contribution of the nuclei energy. This has 
led to the 'static', topological analysis of chemical reactions, and also to an endless debate 
between two communities: PES explorers and dynamicists. The second scheme is based on 
quantum-classical correspondence, i.e., on comparing the differences displayed by an observ- 
able when computed classically and quantum-mechanically. Any deviation between these 
two calculations is then regarded as a signature of 'quantumness', this being a standard 
criterion to discern whether the observable is affected or not by quantum effects. This argu- 
ment, however, can be misleading: the classical and quantum values could be quite similar, 
but the dynamics leading to them could be very different. Therefore, analyzing quantum 
effects in terms of individual events (trajectories) would be highly desirable, which can be 
achieved within a theoretical framework based on the Hamilton- Jacobi (HJ) formalism. 

In the analysis of PESs, molecular processes are described by means of reaction paths 
(RPs), which are continuous curves on the PES connecting two minima through a saddle 
point (first-order maximum) and usually associated with steepest- descent curves The 
minima are usually related to the reactants and products states, while the saddle point 
describes the reaction transition state. All these elements provide us with a mechanistic 
description where the chemical reaction consists of a sequence of RP points, each one related 
to a particular molecular geometry. However, molecular geometry changes as the reaction 
evolves in time cannot be explained satisfactorily, what has derived in dynamical approaches 
and, therefore, the first type of comparative scheme. In this regard, former connections 
between the RP formalism and the classical HJ equation can be found in [sl, for instance. 



The second comparison scheme arises when we realize that the fundamental parameters 
characterizing chemical systems (masses, energies, density distributions, etc.) very often 
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lie on the border between purely classical and quantum behaviors. Actually, many times, 
classifying such behaviors results a hard task and one has to resort to previous experience 
before labelling them as classical, quantum or 'in between'. In this sense, both formal 
0. B Q a.d computafoua. analyses ca. be found ftpcuent.y „> the hterature. the latter 
usually considered simultaneously with the former P, |8[ . The problem is that the solutions 
obtained from classical and quantum evolution equations are formally and conceptually dif- 
ferent. The classical output is a trajectory, a time-ordered list of positions and momenta; 
the quantum output, however, is a wave function and, therefore, a time-ordered list of prob- 
abilities covering the whole position or momentum space. Thus, before comparing classical 
and quantum outputs, some averaging of the former is required in order to approach the 
statistical nature of the latter (semiclassical mechanics constitutes a mixed approach, where 
a classical-like view of quantum effects is achieved after some quantization scheme is carried 
out). The calculation of reaction probabilities and/or cross sections in molecular collisions, 
or Franck-Condon factors in molecular spectra are well-known examples. Though classical- 
trajectory-based studies (and methodologies) are very common within the molecular reaction 
literature, similar quantum trajectory treatments can also be found but scarcely. For in- 
stance, Muga et al. j^, [lo| have computed quantum moments from an averaging procedure 
and then compared with classical quantities. Within the framework of Quantum Hydrody- 



namics 



i.e., treating the probability density as a quantum fluid, chemical reactions were 



brmer y studied by McCuUough and Wyatt [12 , 



13 



14l |. Similarly, in Bohmian mechanics 



ISl . Ilq . the evolution of the system is described in terms of quantum trajectories, which 



can be obtained through any of two approaches |17|: synthetic [17|, Il9| and analytic 
{20! . I21I . [2^. Within the former, quantum trajectories are computed after solving simulta- 
neously the quantum HJ equation and the continuity one, while in the latter (used in this 
work) one starts with Schrodinger's equation and then the trajectories are obtained from 
the phase of the wave function. 

In this work, we explore and discuss the connection between topology and dynamics in 
chemical reactions within a generalized HJ framework. This general scenario is introduced 
in Sec. [TTl Though the corresponding equations are conceptually different, their underlying 
common formal structure allows us to establish a connection among them. This is illustrated 
in Sec. Illlt where we study the dynamics associated with the Miiller-Brown PES model 
23|, I2J], which can be used to simulate, for example, the passage from reactants to products 
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in isomerization or enzymatic (e.g., Michaelis-Menten) reactions [25|]. Classical and quantum 
reaction probabilities are calculated from the corresponding trajectory statistics. Then, a 
comparative analysis between classical and quantum trajectories is presented, as well as the 



26|. Finally, the 



relationship between the latter and the RP through Caratheodory plots 
main conclusions from this work are summarized in Sec. IIVI 



II. A COMMON THEORETICAL FRAMEWORK TO ANALYZE CHEMICAL 
REACTIONS 

The question of whether there exists a common theoretical ground to the geometric 
(topological), classical and quantum descriptions of chemical reactions emerges naturally 
when comparing the steps (1) and (2) above. A priori, they imply two distant mathematical 
frameworks [1]: the PES description involves analyses strictly based on topology, while 
the dynamical ones (classical or quantum-mechanical) are grounded on the HJ formalism. 
However, a closer inspection shows that both frameworks are not so far one from another. 



As shown in the literature [27], within the so-called reaction path models 28|, |29|], PESs can 
be described by an expression similar to the classical HJ equation. Now, based on the fact 
that both the PES topology and classical dynamics are endowed with a HJ equation, it was 



suggested long time ago jSOj the possibility to associate a wave propagation within the HJ 
formalism, this establishing a link with quantum mechanics. Therefore, it should be possible 
to understand the three frameworks on similar grounds, in particular, through the so-called 



Huyghens construction [30^, very well known in Optics [3l|], for instance. This construction 
is a geometric method for mapping the progress in terms of some parameter of surfaces of 
equal phase, which satisfy at each point a HJ-like equation. RPs in the topological case 
and trajectories in the dynamical one then constitute the set of solutions (characteristics) 
crossing transversally a family of such surfaces with the same phase (or equidistant). 

To understand the analogy between RPs and trajectories, note that both can be obtained 
from similar variational principles. That is, consider the functional 

J[x(r)]= r Fi^,±,T')dT', (1) 

J TO 

where x = {xi,X2, ■ ■ ■ ,Xn) denotes a coordinate vector and x = d'^^/dr its derivative with 
respect to an independent parameter r (physically, this parameter is related to the reaction 
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coordinate, s, in the RP formalism and to the time, t, in the dynamical ones). RPs and 
trajectories are curves x(r) for which the integral ([1]) becomes an extremum (maximum or 
minimum) under small variations of the coordinates (i.e., 61 = 0). This gives rise to the 
well-known Euler-Lagrange equations, 

dF _ d dF 

dxj dr dxj 



0, 



1,2 n, 



(2) 



which are satisfied when F(x, x, r) is evaluated along the solution x(r). Alternatively, by 
means of Legendre transformations, it can also be shown [30] that any of these solutions is 
also a solution of a HJ-like equation. 



81 1 

G(x, Vx/, dl/dr, r) = — + -Vll- V^I + V(x, r) 



0, 



(3) 



dr ' 2 

where V(x, t) is a general potential function. The intrinsic structure of -F(x, x, r) is different 
within each one of the three frameworks. However, some similar elements (characteristic 
equations) can still be found among them, as seen below. 

In the RP formalism, the associated F is a homogeneous functional of degree one with 
respect to x j^, i.e.. 



F(x,x, s) 



(x)] 



1/2 r^T 



11/2 



(4) 



where g(x) is the PES gradient at x [32| and the superscript T denotes the transpose vector. 



The corresponding HJ-like equation ([3]) can then be expressed 

(x) ■ g(x) -- 



^ V^r(x) ■ V.V(x) ^"^^ 



3|] as 
0, 



(5) 



where the term 81/ ds is lacking because the PES (here, playing the role of the surface /) 
does not depend explicitly on s, but only on x. From ([5]), 

g(x) = Vxl^(x), (6) 

which, after integration over s, gives the steepest-descent curve, x(s), describing the RP, 
the curve joining two PES minima through a saddle point. This equation is analogous 
to that found in Optics to determine the optical path followed by light in media with 
variable refraction indexes, according to Fermat's principle 3l||, or the geodesic equation in 
Gravitation 33 1. 

In classical dynamics, Fermat's principle translates into Hamilton's principle, and the 
integrand of ([1]) is just the Lagrangian function [3J], which reads as 

1 



F(x,i,t) = ^i^-i-r(x) 



(7) 



in reduced length units (i.e., ^/ ^/m — > x). The corresponding HJ equation ([3]) is 

^^^^^ + \ V^5el(x, t) . Vx5e,(x, t) + F(X) = 0, (8) 

where its solution, the classical action S'd, follows the Huyghens construction [0. Usually, 
in the kind of process that we are interested in, the total energy E conserves and, therefore, 
—dSci/dt = E. This allows to reexpress (IHl) as 

^ V^^ci(x, t) ■ V.5ei(x, t)-[E- V(x)] = 0, (9) 

which has the same formal structure as Eq. ([5]). Similarly to RPs, classical trajectories are 
now obtained from 

i = p(t) = Vx5'ci(x,t) (10) 

by integrating over time with initial conditions Xq and pg. Note that, from this relation and 
Eq. (Q, we can readily derive the well-known Newtonian expression for the velocity 



2[E-V{^)]. (11) 

Unlike the steepest-descent curve, given a PES there is an infinite number of associated 
classical trajectories, as many as initial po one can provide for a given xq and E. Thus, 
in chemical reactivity problems, a single trajectory results meaningless by itself; in order 
to extract valuable information about the process, a distribution of them (i.e., a sampling 
over initial conditions) has to be considered. This statistical problem, equivalent to consider 
an ensemble of identical, non-interacting systems described initially by some pre-assigned 
density distribution, pci(x, 0), can be expressed in terms of a Lagrangian density 13] as 



F{Pch VxPci, Pel, 5'ci, Vx5'ci, S'ci; x, t) 



^ V'^^ci(x,t) ■ VS'ci(x,t) - r(x) I pci(x,t), (12) 



where the variables are now the fields pd and Sd-, and the independent parameters are x and 
t. Applying the Euler-Lagrange equation to (fT2|) with respect to the field variables yields 
(El) and 

+ ■ [pc,(x, t) VxS(x, t)] = 0, (13) 

which is the classical Liouville equation describing a swarm of single, non-interacting parti- 
cles each one evolving according to ([8]). From (fT3|) . note that the time evolution of trajecto- 
ries is totally independent of the density distribution evolution. This can be inferred from 



Eq. (fT2|) . which is separable in Sc\ and pci- The only relationship between pd and individual 
trajectories is that the former is used to choose initial conditions for the latter, this being a 
statistical rather than a dynamical relation. 

It is known |l6, 3^ that the wave formulation of quantum mechanics can be derived from 



a Lagrangian density 

F(^,Vx^,^,^*,Vx^*,^*;x,t) = y 



dt 



dt 



- ^ V^*(x, t) ■ Vx^*(x, t) - l^(x)*(x, t)vl/*(x, t) 

(14) 

when we require the corresponding integral to be stationary with respect to variations in 
the complex- valued field variables and By applying the Euler-Lagrange equation with 
respect to we obtain the time- dependent Schrodinger equation, 



ih 



V2vI/(x,t) + \/(x)vI/(x,t) 



(15) 



dt 2 

(and its conjugate complex if \E'* is considered), where = ■ Vx- However, there is 
an alternative way to proceed consisting of switching from and \E'* to the real-valued 
field variables 5* and p according to the transformation 15|, l35|: \E'(x, t) = p^/^(x, t)e*'^^^'*'*/'^. 
Substituting this expression into (IT^ readily yields the new Lagrangian density [35|, 

'1 



F{p, VxP, p, S, VxS, 5; X, t) 



2V;^5(x,t)- Vx5(x,t)-\/(x) 
/i2v^p(x,t)-VxP(x,t) 



p(x,t) 



8 p(x,t) 

Now, after inserting F into the Euler-Lagrange equation, we obtain 

+ iv^^(x,t)-Vx^(x,t) + Kfr(x,t) = 0, 



(16) 



^P(x, t) 
dt 



+ V^[p(x,t)Vx5(x,t)] = 0, 



(17a) 
(17b) 



where the effective potential 



Kfr(x,t) = \/(x)-^ 



lVjp(x,t) Vxp(x,t) V2p(x,i) 



(18) 



2 p(x,t) p(x,t) p(x,t) 

is a sum of the PES and the so-called quantum potential [l^. Equations fll7al) and (b) 
are, respectively, the quantum HJ equation and the quantum Liouville equation, which are 



coupled, this being the origin of quantum effects. Note that, unhke their classical analogs, 
quantum trajectories undergo a (nonlocal) dependence on their distribution: though fll7bp 
describes the evolution of an ensemble of non-interacting particles (there is no physical 
potential, like V, acting among them), the quantum potential mediates a sort of information 
exchange among them which strongly determines their dynamical evolution. Similarly to 
the classical case, quantum trajectories are obtained from 

x = p(t) = Vx5(x,t), (19) 

with initial conditions Xg, ^(x, 0) and p(x, 0); here, the initial momentum is predetermined 
by the initial phase of the wave function and, therefore, unlike classical trajectories, we 
cannot choose it arbitrarily. Although S" is a multivalued function (i.e., S"(x, t) = S{'x,t) + 
27mh, with n being an integer), this does not affect the calculation of trajectories, since only 
its gradient is needed. Only when p vanishes, this property plays a fundamental role, for it 



rules the appearance of vorticality [16|, l20|, l21| . Moreover, as before, 5'(x, t) and its gradient 



vector field, Vx>S'(x, t), constitute the basic elements for the Huyghens construction. 

The generalized formulation discussed here is summarized in Fig. [1] As seen, math- 
ematically RPs and classical/quantum trajectories are both characteristics coming from 
similar differential equations. However, important differences arise when the terms of the 
corresponding equations are identified. Thus, for a given PES, the RP follows a unique 
geodesic, steepest-descent curve, while there is an infinite set of classical/quantum trajecto- 
ries solutions of the same differential equation (ITOl) / (fT9|) . each one for an initial condition. 
Nonetheless, there is a common point: if one computes a set of steepest-descent curves for 
some given initial conditions xq and go, one fully 'builds' the PES. Similarly, from the set 
of classical trajectories for some given initial conditions E and xq, one will then 'build' the 
action hypersurface 5'c;(x, t) (the same for quantum trajectories). Formally, all the charac- 
teristics are obtained as solutions of ([3]), i.e., Eqs. IQ, (fTOj) and (fT9|) . which present similar 
functional forms. However, while r represents the RP reaction coordinate and for each r 
we have different values of the PES, it plays the role of time for a dynamical (classical or 
quantum) trajectory, which means that at every time we have a different value of the action 
hypersurface. 
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FIG. 1: A schematic view of the generahzed HJ framework presented in this work. 
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III. NUMERICAL SIMULATIONS 



A. The model 

As a working model, we consider the Miiller-Brown PES [23"], which describes the 
reactants-to-products passage through an intermediate (pre-equihbrium) state, 

R^I ^P. (20) 

This PES has three minima, Mi(-0. 558, 1.442), M2(-0.050, 0.467) and M3(0.623, 0.028) 
(within this work, all magnitudes are given in atomic units), which describe the prod- 
ucts, intermediate and reactants states, respectively. It has also two transition states, 
TS'i(— 0.822, 0.624) and TS'2(0. 212, 0.293), which separate products from pre-equilibrium 
and the latter from reactants, respectively. All these energies are indicated in Fig. [2|i along 
the RP, which is described in terms of the arc-length 

N N 

s{x,y) ^ J2 V^^? + = v^(^^ " ^^-i)' + - ^/-i)'' (^i) 

1=1 i=l 

where (xo,?/o) = {xN,yN) = {x,y), and the final point of the RP is Mi. Red circles 
represent the energies associated with minima {Vmi = —0.147, = —0.081 and Vm^ = 
—0.108) and blue ones with maxima {Vrsi — —0.041 and Vts2 — —0.072). A contour-plot 
of the Miiller-Brown PES with the RP (green line) is displayed in the inset. 
Regarding the dynamics, we have considered an initial Gaussian wave packet, 

where Aq = {2Tra,j.ay)~^^'^ , with o"^ = cr^ = ctq = 0.0125. This wave packet describes a 
proton transfer process and, hence, m = 1,836 in the Schrodinger and Bohmian motion 
equations. In all calculations, the initial position of the wave-packet center is kept fixed at 
{xo,yo) = M3, and only its initial momentum, chosen in all cases as {px,o,Py,o) = {—Po,Po), 
has been varied (indeed, the value of po)- The initial conditions for the quantum trajectories 
are obtained randomly by sampling po = |^oP^(Px ~ Px,o)S{Py — Py,o), while for classical 
trajectories we consider this distribution as well as the Wigner one associated with (!22l) . 

Pw{x,p.; y,Py) oc e-(--«)'/2-S-^(p--p-.o)V?i^-(j/-j/o)V2<x?-<x2(p,-p,,o)V/i2_ (23) 
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a) 




Po 

FIG. 2: (a) Transition state energies {TS; blue circles) and reactants, pre-equilibrium and products 
energies (M; red circles) along the RP for the Miiller-Brown PES. In the inset, a contour-plot of 
this PES with the RP (green line); black/red contours represent positive/negative equipotential 
contours, (b) Energy diagram as a function of the initial momentum (see text for details). 

Fig. [2b shows an energy diagram as a function of po, useful to understand the role played 
by the choice of initial momenta in classical/quantum trajectory simulations and, therefore, 
in the calculation of reaction probabilities (see Sec. IIIIB|) . In this diagram, the energy 
of the transition states and minima are indicated by color horizontal lines (blue and red, 
respectively). The black solid line represents the quantum expectation or average value of 
the energy, E = (H), which coincides with the classical average energy when the (classical) 
initial positions and momenta are chosen according to ( l23i) . since 

E= E{x,p^;y,py)pw{x,p^;y,py)dxdydp^dpy = + V + 6 ^ — ^E,{xl,yl), (24) 

where V = YliLi ^ii^o^ Uo)/^ ■> ^ — fi'^/4:ma'^ is the spreading ratio js^, and the last sum runs 
over all (classical/quantum) particles considered (with {xQ,yo) denoting their corresponding 
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initial positions). The black dashed line is i^^, but without S, which coincides with a 
classical ensemble distributed initially according to po- In classical mechanics, the spreading 
ratio can be, therefore, related to the initial momentum distribution around {Px,OyPy,o)- 
Finally, the black dash-dotted line gives the energy of a classical particle initially located at 
M3 and with initial momentum {—po,Po)- 



B. Trajectory ensemble analysis of reaction probabilities 

The analysis of reaction probabilities can be carried out by defining the reaction proba- 



37| 



bility in terms of the restricted norm 

V{t)= / \'^{x,y,t)\'^dxdy, (25) 

where S denotes the space region above the line yR^p{x) = 0.8024a; + 1.2734, here chosen 
as the frontier line separating products from the pre-equilibrium/reactants region (more 
refined choices could be considered, but no significant discrepancies are expected regarding 
products formation at this level). Apart from the amount of products in time, the initial 
slope of this function also provides information about the reaction velocity. 

From a Bohmian viewpoint, fl25p reads as the number of trajectories Nj^ penetratirig into 



E at a time t with respect to the total number of trajectories initially considered 37| 



Wit) ^ (26) 

which approaches V{t) when N ^ 00 and the initial conditions are sampled according to 
Po- Classically, fl26l) can also be applied, but Wci{t) meaning a classical products fraction. 



Probability can flow backwards from products to reactants [12|, which may lead to a multiple 
crossing of ?/j?_+p(x) by the same quantum/classical trajectory. Hence, another interesting 
quantity is the fraction of trajectories going from reactants to products without taking into 
account their return to reactants, 

mt) - (27) 

which, a.t t 00, provides the total amount of products for a given initial state. That is, 
assuming that one could extract the products formed during the reaction by some chemical 
or physical procedure, VV(t) would provide the maximum amount of products. Note that 
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FIG. 3: Reaction probabilities W (solid curve) and W (dashed curve) for: (a) po = 4 and (b) 
Po = 10, and three different initial distributions: quantum trajectories distributed according to po 
(black) and classical trajectories distributed according to po (red) and pw (blue). In panel (c), VV 
at t ^ 700 as a function of po^ black circles denote VVBohm; I'ed triangles, VVg^j^j^, and blue inverted 
triangles, W^ig^er- 

this information is only available when working with trajectories, since one can visualize 
each individual process and, therefore, detect when one reactive event has taken place. 

In Figs. El^a) and[31^b), we have plotted W (solid line) and VV (dashed line) for pq = 4 
and Po = 10, respectively, and three different initial distributions: quantum trajectories 
distributed according to po (black) and classical trajectories distributed according to po 
(red) and pw (blue). Results from (125|) have not been plotted since they are identical to 
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those obtained from (126|) . In all total of 5 X 10^ trajectories has been considered 

and a propagation up to t = 700, the asymptotic value at which VV stabilizes. As seen in 
Fig. [Jb, for po = 4, E — 6 (dashed curve) is well below TS2 and TSi. Therefore, according 
to standard quantum mechanics, the actual dynamics in Fig. [3^ should then proceed mainly 
via tunneling. We clearly observe that WBohm is below W^ig^ier closer to WBohm- This 
can be explained as follows. The classical distributions explore many initial conditions, 
which eventually may imply individual trajectory energies higher than the transitions states 
and therefore lead to the formation of products. This effect is more striking in the case of 
the Wigner distribution than in the classical Bohmian one, where the initial momentum is 
fixed. On the other hand, for Bohmian trajectories the dynamics is radically different: for 
low po, the wave packet spreads faster than it moves and, hence, the expected formation of 
ripples by interference is going to hinder the passage of Bohmian trajectories to products. 
Note that in Bohmian mechanics tunneling does not exist as it is commonly understood in 
standard quantum mechanics, but as an effect arising from the modification in time of the 
barrier due to the quantum potential; a time-dependent, effective barrier, given by (|T8i) . is 
what actually rules the trajectory dynamics. For po = 10, however, E — 6 is above TSi 
in Fig. [2}d, which suggests (quantum-mechanically) a larger amount of products by direct 
passage rather than tunneling. This will cause that the classical Wigner distribution and the 
(quantum) Bohmian one will render much closer results, as it is observed in Fig. [8)3, where 
VyBohm approaches W^igner; ^ut is far from WBohm- this case, though tunneling can still 
be active, the direct passage is going to control the dynamics in both cases, classical and 
quantum-mechanical. Observe that, in the second case, the translational motion is faster 
than the spreading of the wave packet and, therefore, more trajectories can be promoted 
to products before interference starts to influence the trajectory dynamics. Regarding VV, 
we find a trend similar to W, but the difference between the asymptotic values of these 
magnitudes increases with pq, which is due to the larger amount of energy available and 
therefore the increase of recrossings (for the classical Bohmian distribution the effect is 
negligible, because of the less energy available after fixing po). 

In Fig. [St, the maximum amount of products, VV, is plotted as a function of po for 
the same distributions as above. As can be noticed, the formation of products is more 
efficient classically than quantum-mechanically until relatively large values of Pq; indeed, 
both classical distributions provide a larger amount of products for po < 4. As noticed in 
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Fig. [2}d, E crosses TS2 at po ~ 3.5; in Fig. [St, after this po value, the Bohmian distribution 
provides larger reaction probabilities than its classical counterpart, the difference between 
them increasing with pq. As Pq is further increased, the Bohmian and Wigner distributions 
start to get closer, though the latter gives rise to a larger probability until pq ~ 8-9, where 
the situation is reversed. Again, this coincides with E crossing a transition state, this time 
the one connecting the pre-equilibrium with products. 

Summarizing, though an ensemble magnitude may have very similar values classically 
and quantum-mechanically, the underlying dynamics can be radically different, which arises 
from both the effects of the quantum potential and the mathematical structure of the corre- 
sponding motion equations (see Sec. [Ill): in Bohmian mechanics the momentum readjusts at 
each instant in order to satisfy the statistical requirements of quantum mechanics, while in 
classical mechanics positions and momenta are not coupled. Next, we analyze in more detail 
the underlying dynamics in terms of common elements within the generalized HJ framework: 
RPs and classical/quantum trajectories. 

C. Quantum-classical trajectory analysis 

As seen in Fig. |2)d, values of Pq between 8 and 10 are critical, since the dynamics is 
developed around the TSi from the pre-equilibrium region. Thus, in Fig. HI two pairs of 
trajectories for po = 9 have been purposely chosen. In each panel, both members of the 
pair start exactly with the same initial conditions. However, in (a) the quantum trajectory 
(green) is reactive through two-dimensional tunneling, while its classical partner (blue) is 
inelastic, whereas in (b) both trajectories are inelastic. The different and richer behaviors 
displayed by quantum trajectories are due to their motion being ruled by the quantum 
potential (apart from the PES), which is lacking in classical dynamics. 

In order to better understand these behaviors, we consider differences AA(t) = Aq(t) — 
Aci{t), where Ag{t) and Adit) denote the quantum magnitude and its classical counterpart, 
respectively. Thus, in Fig. [5l we show the differences in positions (Ax and Ay) and momenta 
(Apa; and Apy) for the pairs of trajectories displayed in Fig. |4] up to the time when one 
member of the pair approaches the first saddle point {t ^ 100 a.u.). Moreover, the total 
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FIG. 4: Classical (blue) and quantum (green) trajectories representative of the dynamics on a 
Miiller-Brown PES (in red, the RP): (a) the quantum trajectory is reactive through two-dimensional 
tunneling, while its classical counterpart is inelastic; (b) both trajectories are inelastic. The classical 
and quantum trajectories are launched from the initial positions indicated at the lower left corner 
of each respective panel and with the same initial momentum, pQ = 9. 

energy, computed from the classical-like expression 

is also plotted along each classical (thinner)/quantum (thicker) trajectory (blue lines) taking 
into account the energy scale in the right vertical axes of both panels of the corresponding 
figure. At t = 0, fl28l) is between Vts2 (see Fig. [2}d) for both classical trajectories 

and, therefore, none of them will be able to overcome the barrier leading to products at any 
subsequent time since, classically, the energy conserves. Hence, they will remain wandering 
within the reactants and pre-equilibrium regions. Regarding the quantum dynamics, the 
trajectory in Fig. basically follows the RP, thus reaching the products region. When 
compared with its classical counterpart, we note that, at the earlier stages of its time evo- 
lution, the separation of both types of trajectories is increasing. This fact suggests that the 
quantum trajectory is visiting a wider region of phase space than the corresponding classical 
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one (termed as dynamical tunneling hj Heller 38|]), even before reaching the potential barrier 
(barrier tunneling). Thus we could say that we have a double contribution to tunneling. The 
main contribution to the reaction probability will come from the dynamical tunneling since, 
in two or higher dimensions, the trajectory can skip the barrier very easily by exploring 
new regions of phase space. On the contrary, the quantum trajectory of Fig. |3)d is inelastic 
because there is no route to products mediated by tunneling. Arguing in topological terms, 
this trajectory is launched pointing towards a part of the PES with a stronger gradient 
(i.e., a more repulsive wall); the trajectory is not capable then to overcome the increasing 
bending and go backwards, so it remains trapped at the pre-equilibrium region. This dy- 
namical picture of tunneling in two dimensions thus seems to be more demanding than in 
one dimension, since not only energy criteria, but also the initial momentum orientation 
should be taken into account. 

Energetically, quantum dynamics look so different from classical ones (for instance, we 
may observe reactivity quantum-mechanically, but not classically, as mentioned above) be- 
cause ( l28l) does not conserve along quantum trajectories. This can be appreciated in Fig. [5l 



when comparing the value of this magnitude for the two quantum 



trajectories along time. 



In Fig. [5^, there is a relatively fast initial expansion or 'boost' 18|, ll9|] which leads to an 
increase of the energy and therefore enhances the possibility of tunneling though the initial 
energy was smaller than the barrier height. On the contrary, the second trajectory does 
not benefit from this expansion (its energy remains nearly conserved) and, therefore, no 
tunneling is expected. 



D. Statics-dynamics analysis 

Finally, we would like to provide some clues aimed to clarify the long-standing static- 
dynamic controversy, i.e., whether statics suffices to establish reaction mechanisms or, on the 
contrary, true reaction mechanisms are very different from the information provided by stat- 
ics and, therefore, dynamics is compulsory. A direct approach to tackle this problem consists 
of studying the differences in configuration space between the RP and classical/quantum tra- 
jectories. In principle, any static-dynamic comparison might be hampered by the different 
physical meaning of the independent parameter characterizing each curve as well as the 
dimensionality of the problem. In order to avoid these difficulties, though at the expense of 
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FIG. 5: Quantum-classical differences for the pairs of trajectories displayed in Figs.Sl^a) and|H^b): 
Ax{t) (black solid), Ay{t) (black dashed), Apx{t) (red solid) and Apy{t) (red dashed). The total 
energy (|28p is also displayed (the scale is on the right side of each panel) , with blue thinner / thicker 
line for classical/quantum trajectories (the classical energy is constant with time). In the aside 
plots, solid lines indicate the portion of the trajectories up to the time when one member of the 
pair approaches the first saddle point (t 100 a.u.). 

losing some dynamical details, one can consider projection schemes. For instance, trying to 



check the accuracy of the RP formalism, Taketsugu and Gordon [39|] tackled such a problem 
using two distances, the perpendicular and the distance of closest approach. We propose 
a generalization of this approach by computing the whole set of distances between the RP 
and a given trajectory. This generates a matrix with elements 

a, = |x(ti) -XRp(s,)p, (29) 

where the (i, j)-element gives the distance between a (classical or quantum) trajectory, x(tj), 
at a time ti, and the RP XRp(sj), at a value Sj of the reaction coordinate (of course, the same 
division in segments is assumed for the trajectory and the RP), and whose visual display 



gives rise to the so-called Caratheodory plot 26|]. These plots have the property that they 



are always three-dimensional surfaces regardless of the number of degrees of freedom of the 
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FIG. 6: In panels (a) and (b), Caratheodory plots associated with the quantum trajectories dis- 
played in Figs. IHa) andS^b), respectively. The transition from red to blue indicates an increasing 
distance between the RP and the respective quantum trajectory. 

problem because they are based on the definition of distance (|29l) . In Fig. El we present 
the Caratheodory plots associated with the quantum trajectories displayed in Figs. Hl^a) 
and 111(b). The quantum trajectory of Fig. is able to tunnel across the highest TS and 
reach the products region keeping its kinetic energy considerably lower than that of classical 
reactive trajectories. Hence, it remains reasonably close to the RP rather than oscillating 
around it (as it might be expected classically), this leading to a Caratheodory plot which 
displays a valley essentially along its diagonal (see Fig. EK). On the contrary, the quantum 
trajectory in Fig. |Dd is inelastic and, therefore, departs remarkably from the RP, this giving 
rise to a minimum valley out of the diagonal (see Fig. (6)3). 
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IV. CONCLUSIONS 



Here, we have shown the guidehnes to carry out analyses of molecular reactions within a 
generalized HJ framework, which allows to encompass both topology and classical/quantum 
dynamics: static RPs arise when a variational kernel is associated with a PES gradient, 
whereas classical/quantum trajectories are obtained when it corresponds to properly de- 
fined classical/quantum Lagrangians. Specifically, the application of this scheme to the 
topology and dynamics in a Miiller-Brown PES has shown outstanding features, such as a 
better understanding of the discrepancies found between classical and quantum statistical 
results. For instance, reaction probabilities have been computed with Bohmian and Wigner 
distributions of initial conditions. As seen, though both rates look pretty close for some 
values of po, their corresponding underlying dynamics are dramatically different because 
of the quantum-mechanical coupling between individual trajectories and their distribution. 
This leads to dynamics totally different from their classical analogs and, therefore, to ef- 
fects not observable classically, such as dynamical tunneling. It is also worth stressing that 
computing VV and displaying its asymptotic value as a function of po (or another parameter 
characterizing the initial state), can be of interest from an experimental viewpoint, for it 
constitutes a simple method to determine the maximum amount of products formed during 
a reaction. Since VVmax provides an upper bound for the reaction probability, it could be 
used to establish some control mechanism on the reaction based on modifying one parame- 
ter in the initial state preparation (here, po, which could be varied by using laser pumping 
techniques, for instance). On the other hand, quantum-classical differences, obtained from a 
one-to-one comparison of classical and quantum trajectories with the same initial conditions, 
give evidence of dynamical tunneling at early time evolution stages, this contribution being 
mainly the responsible for the occurrence of reaction at certain initial conditions. Barrier 
tunneling is expected to be very small. Finally, the construction of Carathcodory plots has 
been shown to be an interesting tool to analyze static-dynamic differences (i.e., between RPs 
and classical/quantum trajectories). Depending on whether the trajectory is reactive or not, 
these plots show distinct topographical shapes: tunneling reactive trajectories lead to diag- 
onally dominated valleys, whereas inelastic, non-reactive ones lead to a much more involved 
shape. Furthermore, since these plots do not depend on the problem dimensionality, they 
could result advantageous to analyze processes involving large molecules. 
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